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1 Introduction 



Diffractive high energy elastic scattering for hadrons appears to find a satisfactory 
explanation in the framework of the Regge theory with the exchange of the Pomeron 
trajectory Rising cross sections can be accounted for by assuming a Pomeron 
intercept slightly higher than one PJ or, in a QCD approach, by considering the 
Pomeron as a gluon ladder ||. The growth of total cross sections can also be 
described in a way compatible with the Froissart bound in the eikonal model || 
or by assuming that the Pomeron is represented by a double pole in the complex 
J-plane [g. 

In contrast with the above picture, inclusive diffractive collisions in proton- 
proton or antiproton-proton scattering, where one of the initial particles changes 
only slightly its momentum and appears in the final state isolated in rapidity, seem 
to require deep modifications to the standard Regge models. The basic problem with 
diffraction dissociation, known for long ||, is that the integrated cross section osd 
appears to grow faster than the total cross section <7t, thus violating unitarity. For 
example, in the case of a supercritical Pomeron with a(0) = 1 + 5, <jsd grows twice 
as fast, ~ s 2S , as the total cross section does, ot ~ s s . The only way to resolve 
this discrepancy seemed to require the vanishing of the triple Pomeron coupling 
(Pomeron decoupling ||), which however contradicts the experimental data 0. 

A number of different unitarization recipes have been proposed in order to mod- 
ify the energy dependence of the predicted cross section. Eikonal corrections |7) 
succeed in reproducing the main features of single diffraction at high energy, while 
the same effect can be reached by the inclusion of cuts in the Regge theory M. Re- 



cently a different, more phenomenological, approach has been considered P, 10, 11 



Renormalization [[5], |TU| or damping [J [| of the Pomeron flux, that consists in setting 
a limit to the probability that the proton emits a Pomeron, allow for a rising of 



the total diffraction cross section compatible with experimental data. All the above 
approaches are based on a supercritical Pomeron input with a Pomeron intercept 
larger than one. 

Apart from the incompatibility with the experimentally rising cross sections, a 
unity intercept Pomeron would present analogous problems with unitarity ||. If, 
however, the partial wave amplitude, for the Pomeron exchange, presents a simple 
and a double pole in the complex J— plane, cross sections will grow with energy and 
it will be possible to satisfy unitarity at the Born level, without eikonalization j|. In 
all (or most) of the models explored until now, the Pomeron was assumed to be an 
isolated single Regge pole. From QCD we know however |2| that the Pomeron is not 
a single pole, but rather a complicated set of singularities in the J— plane. A simple 
and feasible way to approximate this complicated structure is to take the sum of a 
simple and a double pole (dipole). The dipole Pomeron is known [ffl to have unique 
properties since it reproduces itself under unitarization and thus one expects that it 
can be used also to resolve the abovementioned problem in diffraction dissociation. 
Obviously, the sum of a simple and double pole - like any combination of Regge 
singularities - looses factorizability, although each term remains factorizable. Since 
Regge pole factorization appears to be in conflict with experimental results || , the 
approach we consider is favoured. 

The dipole Pomeron model has been tested successfully in elastic hadron-hadron 



and 7— hadron reactions [HI O, [T3, O, |15| and an application to single diffractive 



dissociation has been considered in Refs. [16|, T7\. It turns out that, in this approach 



the Pomeron contribution consists of two terms, one increasing like the logarithm 
of the energy and the other being energy independent, multiplied by "a priori" 
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different t— dependent vertex functions. This feature, and the assumption that the 
Pomeron couples in a different way to Pomerons and hadrons, opens the way to a 
unified treatment of elastic and production amplitudes. 

Since the inclusive process of hadron diffraction has been discussed extensively 
in the literature [|5|, [T8|, |2D|, [H], we will start, in Section 2, from the Mueller 



discontinuity formula and adapt it to the chosen model. The triple Pomeron con- 
tribution will be discussed in detail and the possibility to satisfy the unitarity con- 
straint will be investigated. While a proof of the proposed solution cannot be given 
in the framework of the Regge theory because the t— dependence of the vertices is 
arbitrary to a great extent, plausibility arguments can be advanced on the basis of 
dynamical models for the Pomeron. Section 3 will be devoted to the inclusion of 
secondary Regge trajectories and the final expression for the cross section will be 
compared with experimental data in Section 4. The conclusion of this work will be 
drawn in Section 5. 



2 The triple Pomeron in diffractive dissociation 



Consider first the process a + b —>■ c + X with the exchange of Regge trajectories 
{i}. From the Mueller discontinuity formula [22 we get 



dp c 16ns x 



(1) 



in the usual Regge pole model. M 2 is the squared mass of the unrevealed state X, 
a>i{t) represents the Regge trajectory exchanged and 

1 ± exp(— inoii{t)) 



sin(7TQ!j(i)) 



is its signature. In the following i — P, f, 7T and ou, where P stands for the Pomeron 
trajectory. 

Consider now the elastic scattering and suppose that, asymptotically, the ab- 
sorptive part in the s-channel, A(s,t), goes like 

A(s,t) ot ^{t) (3 2 {t)s a{t) [h{t) In s + C] , 

then the partial wave amplitude presents a simple and double pole in the complex 
J— plane. The amplitude for the Pomeron exchange can then be written as 



'_ is Wt) 



T(s,t)<x- .\ > mm 

sm(7ra(t)/2) 



h(t) (ins +C 



(2) 



where constant terms have been collected in C. The explicit form of h(t) depends on 
the model. As an example, in a dual model, if the residue of the simple pole has the 
form P(a(t)), the residue of the double pole will be given by / (3(a) da + const ||. 
The form of this residue is such that the coefficient of the double pole can vanish 
for t = 0, if this is required from general principles. 



To substantiate this possibility, we can generalize the picture of Ref. |£| and 
suppose that the Pomeron pole couples to quarks through the exchange of two glu- 
ons. In order to describe the gluon-Pomeron-gluon vertex we use, at high energy, 
the rules of covariant Reggeization [24|, |25| for the coupling of the Pomeron pole to 



massive vector mesons, since the gluons are off-shell. This require the introduction 
of five unknown functions of t (some of them will vanish because of gauge invari- 
ance) that will appear in the amplitude multiplied by polynomials in a(t) and by 

appropriate powers of u, v a ^~ n , where v w s/2 at high energy and n — 0, 1, 

The final expression is complicated, but it can be easily seen that there will appear, 
among the leading contributions u a ^\ at least one term vanishing when a(t) = 1. 



This property remains true when we integrate over the gluon and quark momenta 
and, by taking the derivative with respect to a(t), we introduce a double pole for 
the Pomeron. Moreover this result does not depend on the choice of the partonic 
wave function of the hadron. The conclusion is that, also in elastic hadron-hadron 
scattering, the presence of a term vanishing with t, together with other terms finite 
at t = 0, is highly probable in the conventional Regge residue. It could well happen 
that, if we consider the triple Pomeron vertex, the unitarity condition will impose 
constraints on the couplings such that h(t) in Eq. (|2]) vanishes at t — 0. 

In the dipole Pomeron approach, Eq. (|l|) becomes 
d 2 a 



dM 2 dt 
1 



167TS 2 



E 

X 



ap(t) T / S 7T 



^Pb^X I 1 j2 



(M 2 ,t) 



+ E/aMt)F^ X (M 2 ,t) (^) 



di{t) 



(3) 



Let us consider now the triple Pomeron contribution to Eq. @, neglecting for the 
moment all the interference terms and replacing the sum over intermediate states 
by a discontinuity in M 2 , 

2a P (t) 



h ®*W + c ) 2 + T h2{t) 



ImT yb (M 2 ,t,a P (t),t bi = 0) 



(4) 



where, according to Eq. 



ImT Pb = a (M 2 H (0) (A + h(0) InM 2 + \' (M 2 ) af ^^ 1 )g(t) 



(5) 



g(t) being the triple Pomeron coupling. A term, decreasing with M 2 , is present 
in Eq. (|5|) since we consider also the secondary / trajectory in P — b scattering. 



Obviously, if h(0) vanishes, the same will be true for h(0). In the following, a.p(t) = 
1 + a't, a' = 0.25 GeV -2 and the standard form for the residue will be assumed: 
C = exp(6t). 

By integrating Eq. (|5]) over t and M 2 we get the Pomeron contribution to the 
single diffractive cross-section, ctsd- We will now show that the constraint <jsd < °t 
for all values of s requires that h(t) oc (— t) 7 with 7 > 1/2. Without changing the 
asymptotic behaviour of the Pomeron-hadron vertex, we can assume thatQ 

h(t) oc f-F^-rV 7 > . (6) 



The proportionality constant in the expression for h(t) is unessential since it can 
be factorized out in by properly rescaling the constant C. 

The proof becomes simpler if, according to experimental findings [^|, [|, we 
consider the triple Pomeron vertex g(t) as constant and neglect ImT Pb in Eq. (Q). 
Then, setting B = 2(6 + a' ln(s j M 2 )) and y = ln(s/M 2 ), the t integral can be easily 
evaluated and reads 



dt e 



-Bt 







t V \ 2 7T 2 / t \ 27 

+ c) + — 



t+U J 4 \t + 1 



r(2 T + l)*(2 7 +l,2;S) (y 2 + ^) +2yCr( 7 +l)*( T + l,2;S) + ^ , (7) 



where \l/(a,c;x) is a confluent hypergeometric function [26|. 



In order to integrate over M 2 we transform to the variable B whose upper limit 
is, asymptotically, proportional to In s. The integral can be evaluated exactly by 
using the elementary relations for the \l/ function and, in the limit B ~ Ins — > 00, 



1 Other choices are possible as well: the function h(t) oc 1 — exp(7i) has a similar behaviour for 
suitable values of 7. 
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the behaviour of ctsd can be inferred from the large variable estimate for ^ [26 



a SD ~ r(2 7 + l)-^— B 2 - 2 ^ + ... + C 2 \nB , (8) 
1-7 

where dots in the l.h.s. stand for terms with a less singular behaviour when s — > oo. 

We note that the singularity for 7 = 1 in Eq. (H) is spurious; the exact result 
does not present singularities for 7 > 0. Since, in the model considered, a? ~ Ins 
and o"5£) < err, from the first term in Eq. @ we must have 2 — 2 7 < 1. Hence, 
the parameter 7, in general, must satisfy the condition 7 > 1/2. This inequality 
is necessary to avoid terms, violating unitarity, that rise faster than Ins. It is 
important to notice that the triple Pomeron contribution does not vanish at t = 
because of the presence of the constant C . 

3 Non-leading contributions and the differential 
cross section 

From now on we select the hadrons participating the process: a and c are antiprotons 
(p) and b is a proton (p). Later, for the evaluation of the total single diffractive 
cross section <jsd, the process a = c = p and b = p will be also taken into account. 



On the basis of historical fits |]T8|, p7| , the to trajectory can be neglected and, 
since the 7r trajectory contributes in a different kinematical region with respect to P 
and /, interference terms between it and P, f are suppressed. Hence, in Eq. (^) the 
sum over i refers only to /, and the 7r contribution will be chosen as in || [27], [29| 



d 2 a 



dM 2 dt 



I n 2 (-+) / s \2a n (t)-l 

WZW^T GsO G {i)a^{M ) , (9) 
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where 

G(t) = , 
w 2.3 -t 

^ pp /(47r) = 14.6 and a n {t) = 0.9 1. 

The / contribution, and its interference with the Pomeron, must now be con- 



sidered. The approximation suggested in ^T] is based on the assumption that 
the / couples to hadrons in just the same way as the Pomeron. This choice avoids 
the proliferation of free parameters and is justified from the consideration that, 
while the / is required by the data |7], ^9|, |31| , its contribution is small, in per- 
centage, and can be approximated. A Pomeron-Pomeron-Reggeon term larger than 
0.15 asD is excluded by high energy data |!2] and is completely ignored in a recent 
analysis [[|. 

Since the model we consider for the Pomeron is different from the conventional, 
supercritical one, we must take care in choosing an appropriate / trajectory. Fits 



with a double Pomeron pole Jl3|, [14], [Lf| require an intercept a/(0) higher than 



the value, usually adopted, about 0.55 |]|. In a recent analysis [|33|, however, 
for the degenerate a 2 /f trajectory, the result a + — 1 = —0.31 ± 0.05 has been 
obtained by refitting all the experimental cross sections considered in For all 
data, with errors added in quadrature, a smaller value for a + has been obtained: 



a + — 1 = —0.34 ± 0.05. An intermediate value, —0.32 has been used in Ref. [ [29 
The coincidence of a/(0) — 1 with —0.32 obtained in the fit of hadronic cross sections 
within different models for the soft Pomeron should not be surprising; in a limited 
energy range a behaviour s e , for e sufficiently small, can be well approximated by a 
term of the form (u + vlns). 



8 



Let a(t) be the difference between the P and / trajectories. If we set 

a(t) = ap(t) — Oif{t) = a(0) — 5t , 

then typical values, adopted in the following, are a(0) ~ 0.34 and 5 ~ 0.65. The / 
contribution 

na(t) \ 



R(s,t) 



k 



cos 



llhit) . I K(l 

— — sin 



-2a(t) 



will appear in the final form of the differential cross section: 



d 2 a 
dtdM 2 



A 

7T 2 



2(b+a' ln(s/M 2 ))t 



+ —h 2 (t) + R(s,t) 



+ 



1 9 2 



-t) 



(1 + l(M 2 ) a fW~ l ) 



47r47rM 2 {t-fi 2 ) 2 r; U 2 / 
where all the constant factors have been collected in A 



a^{M 2 ) 



(10) 



(11) 



In Ref. p0| a value near 7.8 is quoted for the parameter k, appearing in R(s, t); 
since, however, the expression ([11]) has been rescaled, k is here a new parameter. As 
far as the other parameters are concerned, b will be fixed from p—p elastic scattering 
(e.g. b = 2.25 GeV" 2 , consistent with the slopes used in [g, 0, [TJ, [TJ, [TJ, |TJ) and 
cr^ p (M 2 ) in the dipole Pomeron model can be written as 



a7(M 2 ) = 0.565 + 2.902 ln(M 2 ) + 44.388(M 2 ) Q ' (0) ~ 1 



(12) 



inspired by the parametrization used in |IJ. Since the form of h(t) is determined 
only near t — 0, it is well possible that the t-dependence of the cross-section should 
be corrected. Hence, a different value of b could be required from the experimental 
data, but this possibility will not be considered in the following. 
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4 Comparison with data 



When comparing the model with experimental data, we find two kinds of prob- 
lems. The first one is related to the experimental definition of single diffraction 
dissociation. The great variety of phenomenological models, adopted by different 
experimental groups in order to extract the published data, makes the test of any 
new model difficult. Moreover, integrated cross-sections do not refer to the same 
intervals of M 2 and t, for different experimental analyses. The second kind of prob- 
lem resides in our parametrization and is strongly related to the first one. The 
integrated cross section cannot be given in compact form and, since the overall 
normalization of the data has an experimental uncertainty of 15 %, it is not an easy 
task to determine the parameter 7 only from the t-dependence of the cross sections 
at different energies. 

While the pion contribution can be fixed as in Section 3, the parameters relative 



to the / trajectory are different with respect to those of Refs. p0| , pl| , since the 
Pomeron contribution differs from the one proposed there. We are left with three 
parameters for the Pomeron and one for the /, plus an overall constant multiplying 
these contributions, while the 7r term has no free parameters. 

From now on, we adopt the standard variable £ = M 2 /s, that represents the 
fraction of the momentum of the proton carried by the Pomeron. Using the ex- 
pression (|ll|) for d 2 a/d^dt in our model, we performed a global fit of the data at 



y/s = U and 20 GeV of E396 @ and at ^ = 546 and 1800 GeV of the CDF 
collaboration [32|. All the data were taken from the compilation of Ref. || and are 
at fixed t = —0.05 GeV 2 . The range of £ for the data of E396 has been limited to 
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0.0160 4- 0.1013; in the case of the CDF data, we have considered £ in the range 
0.0064 -=- 0.109 for the data at y/s = 546 GeV and in the range 0.0033 0.0918 
for those at y/s = 1800 GeV. We have found that our proposed model nicely fits 
all the data for a large range of values of the parameter 7 larger than 1/2. This 
weak dependence on the value of the parameter 7 was not unexpected, since the fit 
was performed at fixed t. In the particular case of 7 = 2 (which will be justified 
in the following) the fit gives for the remaining parameters the following values: 
C = 0.9802, A = 1.9080, k = 0.9839 and / = 2.3987, with x 2 /d.o.f. « 0.9. In Fig. 1 
we compare the curve resulting from the fit with 7 = 2 with the experimental data. 
We can see that our model succeeds in reproducing the experimental data at dif- 
ferent values of s. We have checked that choosing a different value for 7 produces 
only little changes of the other parameters, but does not affect in a sizeable way 
the shape of the fitting curves. 

We have then fixed the parameters in the expression for d 2 a/d^dt in our model 
according to the result of the fit at t = —0.05 GeV 2 and have checked how it 
reproduces other sets of data, obtained at different t values. We have considered 
the data of Ref. |5J at t = -0.015 GeV 2 and those of the UA8 collaboration |p] 
at the relatively large value of t = —0.95 GeV 2 . In both cases our curves roughly 
reproduce the data (see Figs. 2 and 3), thus indicating that also the t— dependence 
in our model is quite reasonable. 

Finally, we have considered the total single diffractive cross section <Jsd, for the 
process p(p) + p — > p(p) + X as a function of y/s. We have compared our model 



with the experimental data of |32| , pq , p7j , |38[ , p9| from the compilation given in 



where some data have been corrected in order to obtain the diffraction cross section 
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for £ < 0.05. In order to make the comparison, we have numerically integrated 
our expression for d 2 a/d^dt, with the parameters determined by the previous fit, in 
the region lA/s < £ < 0.05 and t < 0. In Fig. 4 we observe that the result of the 
integration, plotted as a function of a/s, is in good agreement with the experimental 
data over all the range of values of s, including the Tevatron energies yfs = 546 
and 1800 GeV. We must stress here that the choice 7 = 2 is essential: values of 7 
lower than 2, but larger than 1/2 in order to satisfy ctsd < &t, would give a too 
fast growth with s, whereas larger values of 7 would cause an undershooting of the 
data at large s. 

5 Conclusions 

In this paper we have considered the proton diffraction dissociation in the dipole 
Pomeron model. In this model the differential cross section d 2 a /d^dt can be written 
in the form given in Eq. (|TTD . From the theoretical point of view, the result in 
Eq. (0) assesses two important properties that seem to be required by the data || . 
First, the exact factorization, typical of the Regge pole model, is lost in the dipole 
Pomeron approach. Second, for t = the Pomeron and pion contributions are 
independent of s and the scaling with M 2 of d 2 a/dM 2 dt\ t= o becomes exact if only 
these terms are considered. Moreover, we remark that this model respects the 
unitarity condition without decoupling of the triple Pomeron vertex. The total 
diffractive cross section rises as ln(lns), i.e. slower than the total p — p cross section 
that, in turn, satisfies the Froissart bound. 

We notice that, in Eq. ([11]), the triple Pomeron coupling and the Pomeron- 
proton cross section are tangled in the multiplicative constant A together with an 
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unknown scale factorized from the function h(t). Hence the fit of the experimental 
data cannot determine the aforesaid quantities but, at any rate, it represents an 
important test of the model. Concerning the comparison with experimental data, 
we have found that this model gives a satisfactory fit to the experimental data for 
d 2 a I 'd£dt with regards both to the £— and t— dependence. Moreover, for a suitable 
choice of the parameter 7, it well reproduces also the data for the total single 
diffr active cross section and allows to predict a value of about 11 mb at the LHC 
energy yfs = 14 TeV. 

We stress that in our model the one-pion contribution, parametrized in Eq. @, 
has been fixed from the beginning, differently from Ref. ||, where a multiplicative 
constant has been considered in front of it as one of the two free parameters to 
be fitted. As for the / contribution, in our model it is well below the limit found 
by CDF [32|. The discrepancies observed at large £ from the data of Ref. |35 



and of UA8 [|29] could arise from an underestimation of the contribution of the ir 



and from neglecting that of the uo. According to Ref. p9 |, the one-pion exchange 
contribution is only a small part of the total non-Pomeron exchange background. 
Also the approximated treatment of the / could be responsible for the disagreement 
at large £. What we need is a more rigorous method for justifying the t-dependence 
of our parametrization and an extensive study of non-leading contributions. We feel 
that a deeper insight in these problems is important for applications of the model 
to other processes. 
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Figure 1: Differential cross sections d 2 a/d^ dt [mb/GeV 2 ] vs £ for several values of s. 



Data are from E396 [ 33] and CDF |32] , compiled in . The solid curves represent 
the model with 7 = 2 and with the remaining parameters determined by the fit. 
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Figure 2: Differential cross sections d 2 a/dC, dt [mb/GeV 2 ] vs £ for several values of 
s and t = —0.015 GeV 2 . The solid curves represent the model with 7 = 2 and with 
the other parameters determined by the fit shown in Fig. 1. Data are from and 
were not included in the fit. 
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Figure 3: Differential cross sections d 2 a/d£dt [mb/GeV 2 ] vs £ for y/s = 630 GeV 
and t = —0.95 GeV 2 . The solid curves represent the model with 7 = 2 and with 
the other parameters determined by the fit shown in Fig. 1. Data are from [213] and 
were not included in the fit. 
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Figure 4: Total single diffraction cross section cjsd vs ^/s, compared with the pre- 
diction of the model. Data are taken from the compilation of Ref. 0. 
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